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Abstract 

We embed the S-MSSM - a singlet extension of the minimal supersymmetric standard model 

'~~^ with an explicit /i-term and supersymmetric mass term - in a gauge-mediated supersymmetry 

o, breaking scheme. We find that by absolving the singlet of its responsibility for solving the /i- 

' problem, we are instead able to solve the little hierarchy problem. Specifically, we find that even 

Qj with a minimal embedding of the S-MSSM into a gauge- mediated scheme, we can easily raise the 

1-S^^ lightest Higgs mass above 114 GeV, while keeping top squarks below the TeV scale, maintaining 

perturbative unification of the gauge couplings, and without tuning the other parameters of the 

^~~^ model. 
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1 Introduction 

While there are several arguments that may lead one to expect supersymmetry (SUSY) as a 
weak-scale extension of the Standard Model, none is more powerful than the cancellation of the 
quadratic divergences that generate the so-called "hierarchy problem" of the Standard Model. 
SUSY avoids the order-by-order fine-tuning that one imposes on the Higgs mass of the Standard 
Model in favor of a doubling of the particle spectrum. This doubling is not without cost, as 
several problems that are resolved naturally within the Standard Model {e.g., conservation of 
baryon and lepton number, suppression of flavor-changing neutral currents) require additional 
symmetries or structure in order to be accommodated within the minimal SUSY Standard 
Model (MSSM). 

In order for SUSY to cleanly solve the hierarchy problem, the new SUSY partners must 
lie close to the weak scale - as the mass scale of the superpartners grows, fine-tunings of 
0(?7z^/MgugY) reappear. At present, this problem is somewhat academic since none of the 
superpartner masses have been measured. But constraints from experiment have begun to force 
the MSSM into regimes in which some tuning is required. The worst of these tunings comes 
from the LEP lower bound on the Standard Model Higgs mass of 114 GeV |T]. This constraint 
is a problem for the MSSM for two reasons. First, in large regions of parameter space, the 
lightest scalar Higgs boson, h^, behaves indistinguishably from the Standard Model Higgs, and 
so for these Higgs bosons the LEP bound applies. Second, within the MSSM, the mass of h^ 
is bounded from above at tree level by the Z^ mass, mz, and so must be lifted above the LEP 
bound by large radiative corrections. 

The MSSM does allow for large radiative corrections to the lightest Higgs mass, but the 
corrections necessary in order to evade the LEP bound require that the top squarks become 
quite heavy and/or that the trilinear stop mixing term (At) take particular values which are 
themselves very large. These large stop masses/stop mixing terms feed back into the mass 
parameters of the Higgs potential, pulling the Higgs vacuum expectation value (vev) up to 
scales above a TeV. Stabilizing the weak scale at v = 174 GeV requires tunings of better than 
1% amongst the parameters of the Higgs potential. This smaller version of the original hierarchy 
problem is aptly named the "little hierarchy problem" and is quite generic within the MSSM [2]. 

In the last decade, a great deal of work has been done on finding extensions to the MSSM 
which naturally evade the LEP bound without requiring new tunings. These extensions include: 
models in which the SUSY Higgs evades experimental discovery and thus can lie below the LEP 
bound [3j; models in which the cut-off scale is explicitly made to lie near the weak scale L4J or in 
which the cut-off is lowered implicitly by imposing strong coupling on the theory [5] ; models in 
which new operators [H E] or symmetries [8] are grafted onto the MSSM; and models in which 
new particles are added to the spectrum of the MSSM. This last category divides neatly into 
two very different ideas. In one class are those models which add matter fields with the goal 
that these fields couple to the Higgs sector and generate new radiative corrections to the Higgs 
mass [HI E]. The second, and most popular, idea is to add fields to the Higgs sector itself, in 
order to generate new quartic terms in the Higgs potential, pushing up the Higgs mass already 
at tree level [TO] . 

The simplest implementation of this last approach is the well-studied Next-to- Minimal SUSY 
Standard Model (NMSSM). The NMSSM extends the spectrum of the MSSM by a singlet super- 
field, S, and then imposes an additional Z^ symmetry on its superpotential. This symmetry, 



under which ah superfields are singly charged, forbids the usual /i-term of the MSSM, but 
replaces it with an interaction among the S and Higgs fields: 

WnMSSM = WVukawa + ^SHuHd + -S , (1) 

where Wyukawa represents the couplings of Hu and Hd to the quark and lepton superfields, and 
A and k. are dimensionless coupling constants. (Here we take the SU(2) contraction HuH^ = 
eijHl^H^^ = H^H^ — H^H^.) When the scalar component of 5 receives a vev, a /i-term 
is dynamically generated. Because the singlet vev is generated by soft SUSY-breaking mass 
terms in the scalar potential, it is naturally of the order of the soft masses and therefore of the 
electroweak scale, solving the "/x problem" of the MSSM. (See Ref. [11] for an extensive review 
of the NMSSM.) However the NMSSM has a number of difficulties and its phenomenologically 
viable parameter space is highly constrained both by experiment and theory, some of which we 
discussed in a previous paper [12]. In particular, it is difficult for the single S'-field to solve 
both the /i-problem and the little hierarchy problem simultaneously, though it can be done with 
careful choices of parameters. 

In our previous paper [12], we presented an extension of the NMSSM in which we absolved 
the singlet of its responsibility for solving the //-problem and found that it then solved the little 
hierarchy problem in a very natural way. We called this model the S-MSSM, since it is the most 
general version of the MSSM with one singlet superfield addecQ The most important element 
in the S-MSSM was the inclusion of explicit masses for the H^Hfi bilinear and the S superfield 
in the superpotential: 

VFsMSSM 3 I^HuHd + -fisS^. (2) 

In contrast to the MSSM, the S-MSSM was shown to raise the mass of the light Higgs 
boson above the LEP bound at tree-level without the need of a heavy stop spectrum. After 
including the full one- loop effective potential and the leading two- loop corrections, we showed 
that the S-MSSM handily produces light Higgs masses above the LEP bound and even as high 
as 140 GeV. In order to generate these large corrections, fig must be a few times heavier than 
the other mass scales in the theory: for smaller fis, the singlet mixes too much into the light 
Higgs, pulling down its mass; for large /i^, the singlet decouples, leaving behind the MSSM and 
its too-low Higgs masses. We also showed that Higgs masses above the LEP bound could be 
obtained "naturally" with stop masses all the way down to 400 GeV. 

But one must be careful discussing naturalness within a completely low-energy effective 
theory, such as the S-MSSM. Parameter choices that appear natural at low energies may turn 
out to be quite unnatural once that low-energy model is embedded in an ultraviolet theory. 
Conversely, choices that appear unnatural at low energies may simply be clues about the struc- 
ture of a perfectly natural ultraviolet completion. In either case, it is difficult to speak of 
naturalness without a more complete theory. In Ref. [12] we used the stop masses as a measure 
of fine-tuning, arguing that if stop masses were significantly smaller than 1 TeV and there were 
no other obvious tunings required to make the model work, then the S-MSSM is in some sense 
natural. In this paper, we will explore whether the parameter choices made in Ref. [12] could 
be generated within a more complete model, one in which the soft SUSY-breaking masses are 
determined by only a few inputs. 



^We did impose i?-parity on the model, so more correctly, it is the most general singlet extension of the 
i?-parity-conserving MSSM. 



In order to do so, we will take a very natural model of SUSY-breaking, namely gauge- 
mediated SUSY breaking (GMSB), and use the relations of GMSB to constrain our parameter 
space. We will show that one still finds large, and apparently natural, regions of parameter 
space in which the lightest Higgs lies above the LEP mass bound. In the process we will see 
that the S-MSSM, by not trying to do as much as the NMSSM, can be made to fit into a GMSB 
scheme in more minimal and natural ways, with a much larger parameter space available in 
which to find phenomenologically viable models. 

Our paper is organized as follows: In Section 2, we will revisit the S-MSSM discussed in 
Ref. |12j . In Section 3 we will review GMSB and the various ways in which it is implemented 
in the NMSSM. We then explain how the S-MSSM is can be embedded in GMSB. In Section 4 
we show the results of a study of the parameter space of the GMSB-embedded S-MSSM, with 
a discussion in the final section. At the end of the paper is an Appendix in which the details 
of the S-MSSM, its potential and Higgs spectrum are laid out for the reader. 

2 The S-MSSM as a Low-Energy Model 

The basic ingredients of the S-MSSM were laid out in Ref. |12j : we will review them here, with 
additional details in the Appendix. 

The most general superpotential one can write for the MSSM with one additional singlet, 
and which preserves -R-parity, is: 

W = W^.,kawa + {^i + \S)HuHd + \^JisS^ + \kS' + ^S. (3) 

This superpotential contains explicit mass terms both for the HuHd pair (the usual /x-term) 
and for the singlet itself (/is); as well as a tadpole term for the singlet and a cubic self-coupling. 

A few comments are now in order. First, in our analysis of the S-MSSM, we will allow ^ and 
US to take arbitrary values at or below the TeV scale, independent of the vev of S, and so we 
give up any attempt at solving the //-problem (or //s-problem) . This is the key trade-off that this 
model makes: give up the solution to the /i-problem in the hope of solving the little hierarchy 
problem. Of course, one could imagine ways of solving the /i-problem with additional fields or 
symmetries, but we are only interested in the effective theory described by our superpotential 
above. 

Second, a bare tadpole term can be removed by field redefinitions, and even though there 
is no symmetry that explicitly forbids it, the non-renormalization theorem will prevent the 
tadpole term from being generated radiatively until SUSY is broken. Once SUSY is broken, 
the tadpole could reappear, but not in a calculable way. If the tadpole coefficient ^ 3> MsusY, 
then the singlet will develop a huge vev and ruin the electroweak symmetry breaking; if the 
S-MSSM is to correctly model physics at the electroweak scale, this cannot be the case and so 
we do not consider it. If ^ <C -/VfsusYj then the tadpole is irrelevant to our discussion, and we 
again ignore it. Only if .^ ~ -A^SUSY can the tadpole be relevant, but since this requires some 
tuning of the model, and since we don't need the tadpole to play any particular role in the 
model, we do not consider it further. 

Third, while the S"^ term is required in the usual NMSSM in order to stabilize the potential 
in the S direction, our potential is stabilized by the explicit mass term, /is- Because the S^ 
term is no longer required, and because its effects tend to be small anyway, we will take k to 



be effectively zero in tfie analysis that follows. One could reanalyze this model for non-zero k, 
but we don't expect that there is much to be gained in this direction. 
Therefore our superpotential is reduced to: 

W^SMSSM = Wvukawa + (/i + XS)HuHd + -flsS^. (4) 

We call the model represented by this superpotential the S-MSSM. Philosophically, the S-MSSM 
is an attempt to barely extend the MSSM; not only is the particle content minimally extended, 
but the structure of the vacuum will be nearly identical to that of the MSSM. In this sense, 
our philosophy is similar to that in Refs. [HITS]. The authors of Ref. [7] do briefly considered a 
model like the S-MSSM, though without examining its implications in any detail. The authors 
of Ref. [13] examine a class of models which overlaps the one being considered here. 

The superpotential in Eq. Q, along with the soft-breaking terms, yields the following 
potential for the neutral Higgs and singlet fields: 

4°^ = {ml^ + \f^ + XS\')\H'j' + {ml^ + \^, + \S\')\H',\' + ml\S\' 
+ l^isS - XH^H^^l^ + {BsS^ - {B^ + \AxS){HlH^^) + h.c) 

+ \{9'+9"mr-\Hi\'f (5) 

where g and g' are the gauge couplings of 5C/(2)iy and C/(l)y. 

Minimizing the potential with respect to -ff^, H^, and S yields three constraints (details can 
be found in the Appendix): 
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where Vg = {S) and Vu^ = {H^^), with t; = (v^ + v^Y^'^ = 174 GeV and 

/"off = /U + Aus, (9) 

^A'.cff = B^ + \vs{^is + Ax). (10) 

There are several interesting limits in this model. The first is the limit in which /Xg is large 
compared to the other masses in the model. In this case, the potential becomes quite steep 
along the singlet direction, forcing the vev of S to lie near the origin while decoupling the singlet 
from the theory, leaving only the MSSM as the effective theory. We could also consider the case 
in which /i^ is small. In this regime, the singlet mixes strongly with the CP-even Higgs bosons, 
decreasing the lightest Higgs boson mass, and thus making it harder to surpass the LEP mass 
bound. However, one may be able to find regions of parameter space in which the light Higgs 
states are effectively hidden from detection at LEP, thus evading the bound in a way similar in 
spirit to the proposal of Dermisek and Gunion [3j , though we will not discuss that possibility 
here. 



The S-MSSM works to maximize the Ught Higgs mass precisely between the two Hmits 
discussed above. If fig is too large, the MSSM predictions are recovered; if fig is too small, the 
mixing drives down the Higgs mass. This can all be seen most easily by expanding the Higgs 
masses of the S-MSSM in powers of l/fJ-s (see the Appendix for full details): 
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where the first term in Eq. (13) represents the masses of the lightest CP-even Higgs bosons, as 



calculated within the MSSM: 
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and A corresponds to the 0(1/^^) correction to the splitting of the scalar Higgs masses: 
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It is obvious from Eqs. ( 11 )-( 13 ) that at very large ^g, S" ~ H2+iA^, both of which have masses 



~ /is, and the remaining Higgs bosons behave as in the MSSM. 

It is also helpful to consider the lightest Higgs mass in the Higgs decoupling limit, in which 
rrij^o is taken large, since this simplifies our expressions considerably, while also maximizing the 
light Higgs mass: 
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Here we have included both the leading and sub- leading terms in the l//is expansion, both of 
which are required to reproduce many of our results. (Our numerical work in the next section 
is done exactly, however.) 

In principle, one can keep increasing A to obtain arbitrarily large Higgs masses. However, 
if we take the apparent gauge coupling unification in the MSSM as a sign that physics must 
remain perturbative up to scales ~ 10^^ GeV, then at the weak scale we must have A ^ 0.7; see 
the Appendix for more details. Choosing values for A > 0.7 is equivalent to placing a new and 
lower cut-off on the theory, an avenue for raising the Higgs mass that is already known to work. 

Using this model with A bounded by perturbativity, and a simplified set of input parameters, 
we showed in Ref. [I2] that the S-MSSM can generate very large Higgs masses (up to 140 GeV) 
without any apparent tunings, apart from those already implicit in the /^-problem. The model 
itself is surprisingly similar to the MSSM, and because we find /i^ > 1 TeV in our regions of 
interest, the low-energy spectrum of the S-MSSM will be that of the MSSM. The only readily 
observable difference will be that the light Higgs mass in the S-MSSM will be substantially 



larger than one would expect given the masses of the stops, which will presumably be measured 
at the LHC. 

In the next section, we will embed the S-MSSM into a model of gauge-mediated SUSY 
breaking, in order to explore its parameter space in a more motivated and complete way. 

3 Gauge Mediation in Models with a Singlet 

Gauge-mediated supersymmetry breaking is a robust scheme for generating the soft masses 
within low-energy models of SUSY, with a very small number of free parameters determin- 
ing the entire spectrum for the model. (See Ref. |14j for a review.) In GMSB, the effects 
of SUSY breaking are mediated through messenger fields {<I>,<I>} interacting with the MSSM 
through gauge interactions, thus providing a natural solution for the SUSY flavor problem. 
The messenger fields transform as a 5 and 5 of SU{5), thus preserving the appearance of gauge 
coupling unification that partially motivates our interest in SUSY. In GMSB theories, one also 
implements a hidden sector whose main effects are parametrized by a chiral superfield/spurion 
X which acquires a non-vanishing vev in both its scalar and F-components: {X) = M + O'^F 
(we will assume F <C M^). In its minimal form, GMSB then couples the X superfield to the 
messenger fields in the superpotential, but not to any of the usual MSSM fields: 

W = Xm. (17) 

The vev of X yields mass terms for the messenger fields: the scalar messengers acquire masses 
Mib vF while the fermionic messengers have mass M, breaking SUSY in the messenger sector. 
SUSY breaking is then mediated to the observable sector by gauge interactions at one- and two- 
loops. Specifically, soft masses for gauginos and sfermions are generated at one- and two-loops 
respectively. At the messenger scale, M, these are given by 

for gauginos, and 

m|(M) = 2nVc/-^f|^f (19) 



i 

for scalars, where n specifies the number of messenger pairs, i indexes the gauge group, and a^ 
and C/ are the gauge coupling and quadratic Casimir associated with each group, i, and chiral 
superfield, /. In GMSB, contributions to the trilinear ^-terms appear first at two loops, while 
comparable mass terms appear at one loop, so these can be simply set to zero at the messenger 
scale. 

So far, no mention has been made of the /x and B^ parameters in GMSB. There are typically 
two approaches one can take here. The first is to derive /i and B^ from the conditions of elec- 
troweak symmetry breaking, putting them into the superpotential and soft-breaking Lagrangian 
by hand. In doing so, one implicitly assumes that the physics responsible for generating ^ and 
B^ can be segregated from the physics which generates the remainder of the soft-breaking 
terms. Alternatively, one can extend the model to produce the fi and/or B^ dynamically, but 
the numerous difficulties that arise in this approach have come to be called the ^/-B^-problem of 



GMSB. For this reason, phenonienological studies of the gauge- mediated MSSM usually favor 
the first approach. 

The problem is quite different in the NMSSM. Here one has no choice but to solve the //- 
and 5^-problems since they arise dynamically from the minimization of the scalar potential: 
H = Xvs and S^ = XAxVg- The literature on embedding the NMSSM into GMSB is quite 
involved, but we now review a few of the leading ideas in order to give the reader a flavor of 
the issues and tunings that arise. 

The minimal GMSB version of the NMSSM is nothing more than the model described above 



in Eqs. ( 17) - ( |19[ ): a single spurion which communicates SUSY-breaking to a set of messengers 
charged under SU{3) x SU{2) x f/(l), who then communicate it to the gauginos and scalars 
at one and two loops respectively. The superpotential is that of the NMSSM plus the minimal 



GMSB superpotential of Eq. (|17|): 

W = XSHuHd + \ nS^ + X$$. (20) 

Unfortunately, the S-scalar, being neutral under the SM gauge groups, does not receive a soft- 
breaking mass, 771^, until 3-loops. Because m^ sets the scale for the S potential, Vg is suppressed, 
too small to generate the /i-term required by the chargino and neutralino experimental mass 
bounds [15]. A closely related problem is that the superpotential above possesses an approx- 
imate Peccei-Quinn/i?-symmetry under which all the MSSM superfields, as well as S carry 
charge +2/3. This symmetry is explicitly broken by ^-terms, and spontaneously broken when 
the singlet gets a vev. But because the ^d-ternis are small within gauge mediation, the pseu- 
doscalar Higgs behaves as an axion and is too light to be phenomenologically viable. Finally, 
this model possesses a Z3 symmetry which is broken spontaneously at the weak scale, seemingly 
generating domain walls in the early universe that would be difficult to inflate away. 

Because the minimal embedding of the NMSSM into gauge mediation is so fraught with 
difficulties, several non-minimal embeddings have been proposed which are phenomenologically 
viable in at least some portion of their parameter space, and which manage to get light Higgs 
masses above the LEP bound, though at some cost. We will quickly review two such proposals. 

The first proposal couples S directly to the messenger fields in the superpotential: 

W = XSHuHd + \HcS^+X^^-T]S^^. (21) 

Here a soft mass-squared for the singlet is generated already at one loop, but requires extra 
insertions of the SUSY-breaking mass scale, producing a mass-squared for S, 

which is too small to produce the required /x-term. Thus the model would seem to fail just 
as the minimal case. However, as pointed out in Ref. [16j, S must share the same symmetry 



properties as the spurion, X, in Eq. (21). Thus there can be tadpoles in the model which are 
cut off by the messenger scale and which could drive S to get a vev. There are two possible 
tadpoles, one in superpotential, W, and one in the potential, V, which scale as: 



^w ~ -^Flog(AVM2) (23) 



where A is some cut off, presumably near M. Of course, the precise values of the tadpole 
contributions depend on physics at very high scales, but if we take these expressions literally, 
one can choose F, M and t? <C 1 so that both tadpoles are of order the weak scale, driving S 
to get a vev of the same size. This model produces Higgs masses above the LEP bound either 
at low tan/3, large A and very low messenger scales, or at large tan/3 with very large M and 
very heavy stops ( ^ 2TeV). In the first case, perturbative unification of the gauge couplings 
at 10^^ GeV is not preserved without some additional physics at intermediate scales |16j . 
One variation on the previous model is to allow for two sets of messenger fields: 

W = XSHuHd + \kS^ + X($i^i + $2^2) - xS^2^i. (25) 

Notice that S couples only to the bilinear $2^1 , a requirement enforced by a Z3 symmetry; 
this same symmetry also forbids a tadpole in S. The non-diagonal couplings to the messengers 
prevents ml from arising until two loops, producing a soft-breaking mass, and thus a /i-term, of 
the right order jTTj . As in the usual GMSB scenarios, 74-terms for the squarks are not generated 
until three loops, and are therefore quite small. 

This model was studied in detail in Ref . [17) , where the authors found two limiting behaviors 
for the light Higgs. In one case the light Higgs was essentially MSSM-like, with a tree-level 
mass below the Z^. Here the usual MSSM radiative corrections can lift the Higgs mass, but 
because this is GMSB, the At-term is nowhere close to the maximal mixing case, and so one 
must resort to stop masses around 2TeV in order to beat the LEP bound. In the second limit, 
corresponding to large A, A^^, and |m^|, the light Higgs mass can be raised above the LEP bound 
by the presence of the singlet, but not by significant amounts. This second case lives right on 
the edge of the allowed parameter space, with couplings that are almost non-perturbative. 

For the S-MSSM, we will accept the philosophy adhered to by most implementations of 
gauge mediation in the MSSM, namely, we will assume the minimal gauge mediation sector, 
with no direct coupling of the singlet to the messengers: 

W={iJi + \S)HuHd + ^^.52 + Xm. (26) 

As noted above, for the NMSSM this does not work, because it generates too small of a soft 
mass for the singlet and therefore too small of a ^-term; in our case, a small ml is no longer a 
problem, since we include an explicit //-term. We will derive the values of /i and B^ necessary 
to produce electroweak symmetry breaking with the appropriate (chosen) value of tan /3. We 
shall take /i^ as an unknown input parameter; and because it plays almost no role whatsoever, 
we will set Bs = 0. Of course, this model has a /i/S^//is-problem, but it is not really any worse 
than the MSSM when embedded in gauge mediation. 

Furthermore, this potential has no symmetries, approximate, discrete or otherwise, that are 
broken at the weak scale. Thus there is no concern about domain walls or light axions. There 
is a concern about generating tadpoles, but because our superpotential is technically natural 
until SUSY is broken, and because the resulting tadpoles are not calculable and depend on 
unknown details of the full, ultraviolet theory, we will assume they play no important role here. 

With these caveats, we will go ahead and use GMSB to calculate the soft breaking param- 
eters in the Lagrangian of the S-MSSM in the next section. 



4 A Gauge-Mediated Embedding 

Our procedure for finding a point in the gauge-mediated S-MSSM parameter space is as follows. 
We begin by choosing the messenger scale, M. In particular, we do not want the contributions 
from supergravity mediation, which scale as F/Mpi to compete with those of gauge mediation, 
so we insist on M < 10^^ GeV. For each point we also choose a value of tan/3, and given 
that value, we set A equal to the largest value consistent with perturbativity all the way up to 
10^^ GeV. (See the Appendix for details on this procedure.) 

Given tan/3, we know the weak-scale values for the Yukawa couplings, and we run them, 
along with A and the gauge couplings, up to the scale M using their respective renormalization 
group equations (RGEs); we work at one loop for soft masses, and two loops for dimensionless 
couplings. At Q = M we choose a value of F and set the soft-breaking masses for the scalars 



and gauginos using the boundary conditions of gauge mediation, Eqs. (18) - (19), with n = 1 
pair of messenger fields. We also set At = Ax = and rn^ = Bg = at M. Given these 
boundary conditions at M, we run all masses and couplings down to the weak scale. Since we 
are not adding any dimensionless couplings to the NMSSM, the relevant RGEs are those of the 
NMSSM with k = ^^ = [18]. The RGEs for ^ and ^s play no role in this analysis. 

At the weak scale, we choose a value for /i. We then minimize the combined Higgs-singlet 
scalar potential using the one-loop effective potential, solving for v, and adjusting the value of F 
until we obtain u = 174 GeV or we discover that no such vacuum is allowed for that /i. Finally, 
from the minimization procedure we also obtain B^ and Vs- We note that the low-energy values 
for ml and Bg are very small and play essentially no role in the analysis. 

The physical spectrum of the model is then calculated from the soft masses, including 
the masses and mixings among the neutralinos, charginos, and stops, all at tree level. In 
the Higgs sector, the minimization is done using the complete one-loop effective potential, 
including the contributions from the singlet sector. The mass of the lightest scalar Higgs is 
then corrected using the leading two-loop contributions. We add these in quadrature to the 
one-loop calculation, using FeynHiggs [19] to obtain the two-loop correction. These (negative) 
corrections are sizable, but since the dominant terms come from gluinos, which are largely 
unaffected by the addition of the singlet, the error entering our analysis from using an MSSM 
two-loop computation should be quite small. 

For any point in the parameter space, we also apply all relevant experimental bounds. The 
most constraining bounds come from searches for neutralinos and charginos at LEP. Because 
the lightest neutralino, x5) turns out to be the next-to-lightest SUSY particle (NLSP) in the 
regions of parameter space which we are studying, it can decay to a gravitino and photon 
either inside or outside the LEP detectors, depending on its lifetime. For decays outside the 
detector, the neutralino is simply missing energy, and the LEP bound of 46 GeV applies. If it 
decays inside the detector, the LEP bound on missing energy plus hard photon demands that 
mo > 96 GeV. The transition from one case to the other occurs at roughly yF ~ 10^ GeV |14j . 

which corresponds roughly to M ~ 10^ GeV; for yF below 10^ GeV, the tighter mass constraint 
applies. The limit on the lightest chargino, 94 GeV for tan/3 < 40 |20j . was usually redundant 
once the neutralino mass bound was applied. 

We will show our results in several ways. We begin by considering two random scans of 
lOK points, one scan each with messenger scale M = 10^*^ GeV and M = 10^'^ GeV, randomly 



varying other inputs within the ranges: 

2 < tan /? < 6 

2 < /uj/i < 5 (27) 

300 GeV <fi< 900 GeV 

Note that we only consider a small range of tan/3. For tan/3 ^ 2, the tree-level MSSM Higgs 
mass is very small and it is very hard to lift its mass above the LEP bound, even with the 
singlet-induced corrections. For tan/3 ^ 6, the singlet corrections disappear, since the leading 
contributions fall as sin 2/3 ~ l/tan/3 for large tan/3. We also choose fi > since it is clear 



from Eq. (16) that this will maximize the light Higgs mass. The upper and lower bounds are, 
roughly, where the LEP bound on the neutralino begin to exclude the parameter space (lower 
bound), and where the stops are pushed well above a TeV (upper bound). Finally, we choose 
^s to be a multiple of /i, so that /ig is always large compared to the other mass scales present. 
(Choosing ^^ > is perfectly general, since we can do so by a redefinition S.) We do this to 
ensure that the actual mixing of the singlet into the light Higgs is small, since any mixing will 
tend to decrease the light Higgs mass, as we showed in Ref. [12j and as we will see again below. 
The points in this parameter space are chosen randomly using a flat distribution between the 
limits given above. 

Figure [T] shows the results of this scan. Here we plot the mass of the lightest Higgs boson, 
calculated using the leading two-loop corrections, against an effective mass parameter, MsusYi 
defined as the geometric mean of the stop masses: 

MsusY=Km,-j'/'. (28) 

The red points have M = 10^^ GeV, the black have M = 10^° GeV. On the plot we also show 
the LEP bound on the Standard Model Higgs of 114 GeV. It is immediately obvious that the 
great majority of the points considered lie above the LEP bound, for either messenger scale. 
In fact, for MsusY — 400 GeV, a large percentage of the points are already consistent with 
all experimental limits, and by MsusY — 500 GeV, more than half the points are consistent. 
Thus we see that the LEP bound can be surpassed quite naturally, without any tuning of our 
parameters, in the gauge-mediated S-MSSM. 

It is also apparent that there is little difference between the two messenger scales. The 
additional running that occurs in the M = 10^^ GeV case has negligible impact on the Higgs 
spectrum. Thus it appears that our results are fairly robust to variation of the messenger 
scale. This is in fact true, until the messenger scale falls below 10^ GeV, when the experimental 
bound on the neutralino becomes significantly stronger due to decay inside the LEP detectors, 
as discussed above. 

As expected, the upper bound on the Higgs mass is lower than we found in Ref. [12j. There 
we treated the S-MSSM as an effective theory, where all parameters could be chosen to maximize 
the Higgs mass, and we found that we could push the Higgs mass above 140 GeV. Here, with 
the additional structure of gauge mediation imposed on the model, we are not able to achieve 
Higgs masses above 124 GeV. Nonetheless, masses above 114 GeV are quite generic and even 
expected within the parameter space we considered. 

One issue that deserves more discussion is the value of fig ■ Throughout the parameter space 
above, we chose fig to be large compared to the other parameters in the Higgs potential. Clearly 
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Figure 1: Scatterplot of the lightest Higgs boson mass versus the effective stop mass, MsusY, for 
10^ randomly selected points with a messenger scale M = 10^*^ GeV (black) and M = 10^^ GeV 
(red), and within the parameter space given in Eq. (27). The two regions overlaps almost 



exactly, except at very large Msusy- The solid black line is the LEP bound of 114 GeV. 



as ^s — )• cxD, the effects of the singlet must decouple. But as discussed above, at low fis, the 
singlet has a substantial mixing into the lightest Higgs, pulling down its mass. In order to 
obtain a large Higgs mass, we must therefore choose fig within some limited range, bounded 
from above and below. This could be a potential source of tuning in our model, and so we plot 
the effect of changing /ig in Figure [2] 

In this figure, we take three values of tan/3 = {2, 2.5, 3}, and plot the maximum mass of the 
lightest Higgs as a function of /is//U, requiring both stops to lie below 1 TeV. (We vary F/M 
and fi to obtain v = 174 GeV and tti^^ = ITeV.) It is evident from Figure 2 that at Hs ^ 2//, 
the mixing of the singlet into the light Higgs causes a rapid decrease in the Higgs mass, pushing 
it below the LEP bound (solid line). However, for /x^ ^ 2fi, there is a wide range which remains 
consistent with LEP. In particular, the decoupling is very slow as fig increases, so that even for 
Us — 10^ one can find points above the LEP bound, at least for some choices of tan/3. Thus 
we do not consider the choice of //s as introducing a new tuning into the model. Instead, the 
LEP bound is simply telling us that fj,s must lie above roughly 2/i but must remain below a few 
TeV, a natural expectation within SUSY models. 

We can see more explicitly the correlation among fis, AfsusY and the Higgs mass in Figure |3J 
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Figure 2: The lightest Higgs mass as a function of /Us//i. For each hne, the messenger scale 
M = 10^^ GeV, and F/M and ^ are chosen so that mr ~ 1 TeV. Each line represents a different 
value of tan (3 as labelled. The lines end at low jjLg when electroweak symmetry breaking fails. 
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Figure 3: A scan of parameter space for tan/3 = 2 and M = 10^^ GeV, varying /i and /x^/// 
within the ranges specified in the text. All points in the figure are consistent with experimental 
bounds, including the bound on the Higgs mass of 114 GeV. The points are color-coded by 
the light Higgs mass calculated from the model parameters, in 2 GeV steps, beginning with 
m/jO < 116 GeV on the outside of the triangle, and increasing to m/jO < 118, 120, and 122 GeV 
as one moves into the triangle and to the right. 
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89 
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4 


350 
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99 


10^^ 


2 


400 


1000 


360, 460 


115 


74 


10^^ 


2 


850 


2125 


710, 960 


123 


85 


10^^ 


2.5 


750 


2060 


730, 960 


124 


94 


10^^ 


3 


700 


2450 


730, 970 


123 


100 


10^^ 


6 


400 


2400 


470, 630 


114 


107 


10^^ 


6 


600 


3600 


720, 940 


118 


111 



Table 1: A sampling of points with Higgs masses above the LEP bound, either chosen to 
maximize rrii^o or to minimize the stop masses. Ah masses are in GeV. 



Here we have chosen a single value of tan/3 = 2 and M = 10^^ GeV, and varied the other 
parameters as before. We plot the stop mass scale MsusY versus ^s, keeping only those points 
consistent with all experimental bounds, including the Higgs mass bound. We then color code 
the points based on the calculated Higgs mass. Notice that the Higgs mass constraint generates 
a triangular region in the figure. For very light stops, only a small range in /x^ is allowed, but 
as the stops become heavier, the range for fig quickly grows to allow larger values of //«, and 
even heavier /i". However, it does not grow to include smaller fig (and in fact Hs — 1 TeV is no 
longer allowed for heavy stops though it is for light stops), an artifact of the GMSB boundary 
conditions, which implicitly correlate /x to the stop masses, requiring larger and larger /u as 
-^SUSY increases. From the color coding one also sees the behavior seen in Figure [2j where the 
Higgs mass first increases with increasing /i^ but then decreases as decoupling sets in. 

In order to provide the reader with more details about the spectra predicted in these models, 
we have chosen a number of models which were found in scans with M = 10^, 10^*^ and 10^^ GeV 
and displayed their inputs parameters in Table [T] We also show the masses of the stops and the 
calculated Higgs mass. The points with the heavier stop mass just below a TeV were chosen 
to maximize the Higgs mass, though at the cost of heavier stops. The points with light stops 
were chosen to minimize the stop masses, though at the cost of lighter Higgs masses. 

For comparison, we also show the lightest Higgs mass as predicted in the MSSM using 
the same low-energy spectrum; specifically, we take all sparticle masses and m^o = m^o and 
calculate rrif^o in the MSSM. Prom the perspective of the LHC, both of these sets of models 
would look identical, except for the greatly increased Higgs mass in the S-MSSM. In particular, 
notice that all the points in the table would be ruled out by the LEP bound in the MSSM 
alone. 

Though it is not obvious from the table, as one changes the messenger scale M, there is 
very little change in the values of the Higgs mass that one obtains for a given stop mass scale 
(see Figure llj). However, this breaks down for M at or below 10^ GeV case, for which we must 
apply the stricter mass constraint on the neutralinos, accounting for the decay of the neutralino 
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inside the LEP detectors. For both of these reasons, models with heavy Higgs and hght stops 
are very hard to generate for a messenger scale at or below 10^ GeV. 

For completeness, we also take one of the points from the table (the row with the lightest 
stops) and show its complete spectrum. Specifically, for M = 10^^, tan/3 = 2, fi = 400 GeV, 
^s = 1 TeV, we find uif^o = 115 GeV, and: 



Mi,2,3 = {60,125,410}GeV 

rriQ -^ j^ = {495, 460, 455} GeV 

^0 = {60, 110, 410, 420, 1010} GeV 


m^^ 2 = {360, 460} GeV 

mi-^ = {230, 130} GeV 

m+ ={150, 595} GeV 


rr>-AO,H^,H± = {510, 515, 515} GeV 


m^oHO = {1010, 1010} GeV 



Notice that there are now five neutralinos, with the inclusion of the "singlino", the fermionic 
component of S. A number of sparticles in this particular spectrum are easily discovered at the 
LHC. However the Higgs bosons, apart from the lightest Higgs, will be a challenge for the LHC 
and may require a next generation lepton collider. The singlet, which is quite heavy, is well 
beyond any hope for discovery at the LHC, and so its presence must be felt indirectly, through 
the mass correction it induces in the lightest Higgs. 

If this spectrum of squarks, sleptons, neutralinos and charginos were to be discovered at 
the LHC and interpreted in terms of the MSSM, one would predict a light Higgs mass of only 
74 GeV, far below the LEP bound. Even if one were to allow for maximal mixing (i.e., large At, 
in contradiction to gauge mediation), one could push the Higgs mass up, but only to 102 GeV. 
Thus, observation of this model at the LHC, with a light Higgs mass above 102 GeV, would 
require physics beyond the MSSM, and the S-MSSM would have to be a strong candidate for 
that new physics. Unfortunately, the smoking gun of the S-MSSM, the singlet, could hide itself 
from direct discovery for many years. 

Of course, the fact that this is a gauge-mediated model does imply some clear experimental 
effects, including correlations among the sparticle masses and possibly the decay of the NLSP 
inside the detector. From a cosmological perspective, the gravitino is the stable dark matter 
candidate, but our singlet sector plays no role in the calculation of the relic abundances. 



Finally, although the results shown above assume one pair of messengers (n = 1 in Eqs. ( 18 ) 



(19)), we repeated our analysis for n = 2,3 and found no significant differences between the 



three cases. Therefore we only show the n = 1 case here. 

5 Conclusions 

The S-MSSM, which we introduced in Ref. fT2] . is simply a generalized version of the NMSSM 
- it is the MSSM with a singlet, but also with an explicit ^-term and mass term for the singlet. 
While most singlet extensions of the MSSM attempt to solve the /i-problem, we are using the 
singlet to solve the little hierarchy problem. In our previous paper we found Higgs masses as 
large as 140 GeV with stops below 1 TeV. As a low-energy effective theory, however, we were 
free to choose our parameters to fit the problem at hand, namely pushing up the Higgs mass. 
While we did find models that are free of stop-induced fine tuning using sets of inputs that did 
not themselves appear to be tuned, it is possible that we were making choices for the parameters 
that could not be realized in a realistic model of SUSY-breaking, and would therefore be less 
natural than they appear. 
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In this paper, we embedded the S-MSSSM into a gauge-mediated SUSY breaking scheme. 
Within this scheme, there are very few free parameters, and so the low-energy parameters are 
highly constrained. Furthermore, because the A-terms are generated at higher loops and are 
thus small, it is difficult to push up the light Higgs mass using the ^j-dependent corrections. 
Thus gauge mediation provides a much more stringent test of the S-MSSM idea. From our 
results here, it appears that the S-MSSM remains viable once embedded in GMSB, producing 
Higgs masses up to about 122 GeV without fine tuning in either the stop sector or anywhere 
else in the model. The choices of parameters that gave Higgs masses above the LEP bound 
are quite generic. To a good approximation, in order to obtain Higgs masses above the LEP 
bound one only requires that tan/3 is low (between 2 and 6), /Ug is several times fi (with ^ > 0), 
and that the messenger scale is not too low. Most model points that fit this description are 
consistent with all direct experimental constraints. And all this is done without resorting to 
non-perturbative values for A, lowering the cut off, or having stop masses in the multi-TeV 
range, as with some competing models. 

Along the way we found that it is not very difficult to embed the S-MSSM into a gauge- 
mediated scheme, which is quite unlike the NMSSM. We do not require sizable rn^ as in the 
NMSSM, and so non-minimal GMSB embeddings are not necessary. While we do not attempt 
to explain the source of the fi, B^ or ^^-terms, this is no worse than most gauge-mediated em- 
beddings of the MSSM. This implies, that like the gauge-mediated MSSM, the gauge-mediated 
S-MSSM is incomplete. Our primary assumption, however, is that any eventual understanding 
of how II and ^s are generated within this model will not disrupt the scalar potential at the 
weak scale. 

One of the additional benefits of the S-MSSM is that it opens up parameter space that was 
previously unavailable in the NMSSM. In particular, the tunings that are required within the 
NMSSM in order to break the electroweak symmetry while producing a ^-term large enough 
to pass chargino mass constraints are absent here. The model presented here works for almost 
any possible value of A\. And although we set k = A^ = for simplicity, they are both free to 
take on a wide range of values, with only the constraint that k ^ 1 so as not to suppress the 
perturbative bound on A. 

In summary, by adding a singlet to the MSSM, but not requiring that it solve the /i-problem, 
we find a phenomenologically viable, very minimal and quite interesting extension of the MSSM. 
In particular, this model, the S-MSSM, seems to naturally solve the little hierarchy problem of 
the MSSM, and this solution is stable to embedding into a gauge-mediated scheme for SUSY 
breaking. Though the model is almost identical to the MSSM from the perspective of the LHC, 
the enhancement to the Higgs mass, even when stops are light, would be the first evidence that 
this model is the correct description of nature at the weak scale. 
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Appendix 

Here we will give some of the details of the scalar potential and its minimization which were 
left out of the main body of the paper. The model itself it defined by its superpotential, given 
inEq. Q: 

W-SMSSM = Wvukawa + (/^ + XS)HuHd + -f^sS^- (A.l) 

To this we add the corresponding soft-breaking terms: 

V;oft = mljHul^ + mljHal^ + ml\S\^ + {B^HuHd + h.c.) + {AxXSHuHa + h.c.) (A.2) 

The potential for the Hu, H^ and S scalar fields is given by the familiar F- and D-terms, as 
well as soft breaking operators. Written in full, the scalar potential is given by 

V = Vf + Vd + 1/soft 

= {ml^ + l/i + \S\^)\Hu\^ + {ml^ + \^l + A5|2)|F,|2 + ml\S\^ 
+\lisS - XHuHdl^ + {{B^ + \AxS){HuHd) + 5,5^ + h.c.) 

+ \{9^ + 5'^)(|i/„p - m^f + \9'\HiH,\\ (A.3) 

Notice that when S gets a vev, it will shift /x by an amount Xvs. By minimizing this potential 
with respect to the neutral components of the three fields, we obtain: 

"ik +^^^ + ^^^s + A^t;2 cos^ /3 + 2\^ivs + ^{o^ + 9'^){vl - vl) 

= (B^ + {Xfis + XAx)vs) cot /3, 
ml^ +fx^ + X\l + X\^ sin2 /3 + 2XfiVs - \{g^ + 9'^){vl - vj) (A.4) 



B^ + {Xfis + XAx)vs)tan /3, 



^2 sin 2/3,, . , , , , 2 _ / 2 , x2 2 , 2 



where {H^^ = Vu = f sin/3, (-ff^) = f (^ = f cos^, {S) = f^ and v'^ + vj^ = v^ . When these are 
added and subtracted, one obtains Eqs. ([g]) - ([s]) in the text. 

The mass matrix for the CP-odd components of the neutral Higgs and S scalar fields is 
symmetric and given by the following terms: 



(A.5) 



Ml^ -- 


= {Bf, + {Xfis + XAx)vs)cot^ 


Ml^ -- 


= Bf, + {Xfis + XAx)vs 


Mh -- 


= {— fj,s + AxX)v cos P 


Mh -- 


= {B^ + {Xfis + XAx)vs) tan 13 


Mh -- 


= {-fis + AxX)vsinP 


Ml -- 


= nl + X'^v'^ + ml-2Bs. 



Because this matrix has a zero eigenvalue corresponding to the longitudinal polarization of the 
Z^, simple closed forms for the eigenvalues are easy to find. However it is more helpful to 
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express the eigenvalues in terms of an expansion in inverse powers of ^s, in which we assume 
Us > lJ;V,Ax: 

"^Ao = — ^ + 2Aa - — ^ + • • • A.6 

m^o = /uf + (2A2i;^ + m^ - 25^) + • • • (A. 7) 

where the elhpses represent higher orders in 1/ns- (All calculations given in this paper were 
done numerically and without use of these expansions.) 

The mass matrix for the charged Higgs is only 2x2, since the singlet has no charged 
component, but it still contains singlet-induced corrections: 

Mfi = {B^ + {Xfis + XAx)vs) cot 13 + {mlr-X'^v^)cos'^ (3 

Mu = 5/. + (A/x, + AAA)z;s + ^(m^-AV)sin2/3 (A.8) 

MI2 = {B^ + {Xfis + A^a)^^.) tan /? + (m^ - X'^v^) sin^ /3 

using rriyy = g^v^ jl. Besides the massless Goldstone mode, the charged eigenstate has mass: 

Notice that we obtain the usual MSSM result that "i.^± = "w^o + "1h^ in the limit fj,s ^ 00, 
and taking A^ to be the usual A^, and so the heavy Higgs doublet decouples in the large m^o 
limit, just as in the MSSM. 

Finally, the mass matrix for the CP-even Higgs bosons is given by the following terms: 

Mfi = {Bf, + {X^is + XAx)vs) cot 13 + m%sm'^l3 

MI2 = -(-B^j + (A/Us + AAA)fs) -?n,|sin/3cos/3 + 2A^7;^sin/3cos/3 

Mi3 = 2A^t;si' sin /3 - (A/is + A^a)^ cos /3 + 2A;UU sin /3 (A. 10) 

^22 = {Bf, + (A/is + XAx)Vs) tan /? + m| cos^ (3 

M23 = 2A fsWcos/3 — (A//s + A^A)i'sin/3 + 2A/it'cos/3 

Mfg = /is + A^v^ + ml + 2Bs 

The full expressions for the eigenvalues of this 3x3 matrix are much too long to include here, 
but if we assume again that the singlet mass is large and expand in inverse powers of /i^ , then 
the eigenvalues take on a simple form: 



X^v^ 
mlo = mlo + {nsm2f3-Ax-A) + --- (A.ll) 



,2 _ _2 

'mssm Hg 



m'lo = m%o + (/isin2/3-^A + A) + --- (A.12) 



'1 -"-"MSSM 



/^S 



mU = fii + {2X^v' + mi - 2Bs) + ■ ■ ■ (A.13) 



2 



where the terms labelled "MSSM" are the scalar Higgs boson masses obtained from Eq. A. 3 



with X = fXs = Bs = Ax = 0. The mass term A corresponds to the 1//Xs correction to the 
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splitting of the scalar Higgs boson masses and it is given by 

,2 ^2 \„„„2no ..1^2 , ^2 



Ax{mi - m-\o)cos-^2(3 - fi{m^.o + mi) sin 2/? 

A= ^ . (A.14) 

(m^o + m|)^ — 4m^o"^| cos^ 2/3 



In the text we also presented a form for the CP-even Higgs masses, Eq. (16), in which we have 
included the next order in 1/fJ.s while taking ttz^o large. 

Among the parameters that enter the scalar potential, a special role is played by A, the 
coupling of the singlet to the usual Higgs doublets. In this analysis, we always take A to be 
maximal while preserving the unification of the gauge couplings at ~ 10^^ GeV, a constraint 
often called perturbative unification. The maximal value of A must be obtained by running the 
renormalization group equations (RGEs) for A, the Yukawa couplings and the gauge couplings 
between the weak and high scales. We did this at two loops for each coupling. 

The RGEs used in the analysis are the usual RGEs of the NMSSM, though with k = [H]. 
Rather than reproduce all the relevant RGEs here, we will refer the reader to that reference. 
However, the RGE for A itself is worth presenting: 

"-^ _ -^ fifo.^ I „.2\ , „.2 , A\2 3^2 o„2 

dt 167r2 



3(y2 + y2) + y2 ^ 4^2 _ 2 _ 3^2 ^^^^g) 



A /207 4 15 4 022 60, 2 692 

+^5?y? - Igh! + %iA' + i6<7|(y? + y!) - Qyhl 

-9{yf + yt) - 3yt - lOA^ - X\9y^ + 9yl + Sy^) 

where gi = -y/S/S^' is the usual GUT-normalized hypercharge coupling, g2 = 9 and t = logQ 
for arbitrary scale Q. 

The leading terms in the RGE for A are the contributions of top Yukawa coupling and A 
itself: 

Both terms tend to increase A in the ultraviolet and therefore the maximum value of A depends 
sensitively on the value of yt = mi/(f sin/3). Thus the maximum value of A in the infrared, 
Amax, depends sensitively on tan/3. This behavior is evident in Figure |4| where we have plotted 
the maximum value of A as a function of tan/3. Here we define numerically the maximum value 
of X{mw) as that value which runs to A = vr at the gauge unification scale. (The exact value 
we choose for A at the unification scale has almost no effect on our results.) 

Looking at Figure |4j we see that for tan/3 ~ 1, perturbativity breaks down before the 
unification scale for any value of A, but for tan/3 ^ 1.5, non-zero values of A are allowed. 
In fact, the maximum value of A rises quickly at low tan/3 until it plateaus at A ~ 0.7 for 
2.5 ^ tan/3 ^ 40. At very high tan/3, the bottom Yukawa suddenly becomes important and 
Amax is again suppressed. Throughout this analysis we take the value of Amax appropriate for 
the value of tan (3 being studied. 

One last note about taking A = Amax in our analysis. While this may appear to be a tuning, 
it is in fact quite natural. The value of Amax is actually an infrared pseudo-fixed point of the 
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Figure 4: The maximum value of A at the weak scale, consistent with perturbativity of all 
couplings up to the gauge unification scale, as a function of tan/3. 



RGEs - if we choose A to be any large, but perturbative, value at the high scale, it will run 
down to a value very close to Amax in the infrared, regardless of the exact value at the high 
scale. This is a similar behavior to what is observed for the top Yukawa, and in this sense it is 
equally natural as a large yj. 
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